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Abstract 

Quantum random variable, distortion operator are introduced based on canonical operators. As the 
lower bound of rate distortion, the entanglement information rate distortion is achieved by Gaussian map 
for Gaussian source. General Gaussian maps are further reduced to unitary transformations and additive 
noises from the physical meaning of distortion. The entanglement information rate distortion function then 
are calculated for one mode Gaussian source. The rate distortion is accessible at zero distortion point. For 
pure state, the rate distortion function is always zero. In contrast to the distortion defined via fidelity, 
our definition of the distortion makes it possible to calculate the entanglement information rate distortion 
function for Gaussian source. 

Distortion operator: Two major parts in classical information theory are channel capacity and rate distortion 
theory. They concern respectively with the reliability and effectiveness of information transmission. In quantum 
information theory, channel capacity has been widely investigated, but little effort has been put into developing 
quantum rate-distortion theory 1 2 . It was proven (J that the quantum rate-distortion function R{D) is lower 
bounded by entanglement information rate-distortion function (D). For a given source R^iD) is defined by 

R'iD)^ min (1) 

S\d{£)<V 

where d is some distortion function, and £ is the channel. The result is proven under the assumption of distortion 
function defined by transmission fidelity. It is linear among different modes. We can extend the distortion 
function to a more general form. The result will also be true if it is linear among different modes. One of the 
useful distortion function is mean square function as used in classical information theory for Gaussian source. 
The mean square distortion in classical theory is d = E{d{Y,Y'j) = Jp{y,y')d{y,y')dydy', with d{y,y') = 
{y " v'Y ■ Where Y is the input random variable and Y' is the output, y') is the joint density distribution 
function. The same idea should be extended to quantum information theory. What is the quantum corresponding 
of random variable? We prefer the canonical operators X and P. Then the distortion operator will be introduced 
as 

1 " 

B) = - Yy,^A. - XB^f + (Pa. + PBifl (2) 

Where A is the sender and B is the receiver. The Schmidt purification of the sender state pA is obtained by 
introducing the reference system, described by Hilbert space H_r, isomorphic to the Hilbert space Jiq — TIa 
of the initial system, Then there exists a purification of the state pA, a unit vector £ Tig ® Ti-B, such 
that pA = Trii; 1 I .Where 1 1/)) = '^i^\f\k\^k)\^k), with and \\k) are eigenvalue and eigenvector of 
PA respectively. After transmission, the joint state will be p^'^ = {£ ®1) |?/') (?/'| , it is easy to verify that 
Pa = TrQ{p^^ ) and pB = Trii[p^'^ ) = £{pa), which means the system R remains at the input state pA and 
system Q evolves to the output state pB- The average distortion will be 

d^Trp^Q' d{A,B). (3) 

A similar quantity was introduced to obtain entanglement of formation of symmetric Gaussian states The 
average distortion possesses some kind of EPR-uncertainty of the joint state, we here neglect the mean of 
canonical operators for simplicity, thus all of the first moments of the states will be neglected in the follows. 
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Entanglement information rate-distortion function: The coherent information Ic{p,£) — S{£{p)) ~ S{p^'^ ). 
The Entanglement information rate-distortion function is the global minimum of the coherent information 
under the certain distortion. There is a useful lemma in classical information theory which gives necessary and 
sufficient conditions for the global minimum of a convex function of probability distributions in terms of the first 
partial derivatives. The lemma was extended to quantum information theory |3| in evaluating the capacities 
of bosonic Gaussian channels. Let F he a convex function on the set of density operators which contains 
Pq and p, the necessary and sufhcient condition for F achieves minimum on po is that the convex function 
F{{l'-t)po + tp) of the real variable t achieves minimum at t = for any p. That is ^ \t=oF{{^ — t)po + tp) > . 
Here Ic{p,£) is a function of density operator p^*^ . Coherent information is convex due to channel operation[I], 
that is for operation £x = X8i + {1 — X)£2, where < A < 1, one has Ic{p,£\) < XIc{p,£i) + (1 — X)Ic{p,£2)- 
Thus coherent information is a convex function of density operator p^^ . Suppose the minimum is achieved at 
RQ' 



Po 



{£q I) ItJj) {^P\ , the necessary and sufficient condition will be 



^l*=0/c(p,(l 



t)£o + t£)>0 



(4) 



The derivative will be— Tr(£(p)— £o(p)) ^og£Q{p)+Tr{p^'^ ^Po^ ) ^ogp^^ . li£o is a trace preserving completely 
positive (CP) Gaussian operation, then for a Gaussian input state p, the output state £o{p) and the joint 
state Pq^ will be Gaussian. Hence their logarithms are quadratic polynomials in the corresponding canonical 
variables|3]. The derivative will be zero under the constrains of the first and second moments. Where the trace 
preserving property of £ is also used. The conclusion is that for any channels with the same first and second 
moments, Gaussian channel achieves the minimum of coherent information. The moments of the channel is 
with respect to a given Gaussian input state. 

Gaussian channel: Gaussian CP maps are defined as maps which transform Gaussian states into Gaussian 
states. Gaussian CP map is thus isomorphic to bipartite Gaussian state[5]El 



G = [ dxexp{-]x^rx + iD'^x - C)W{x) 



(5) 



where W{x) = exp[—ix'^R] are Weyl operators and R — {Xi,Pi,X2, ■ ■ ■ ,P2n), with [Xfc,P;] ~ iSki-We in the 
following will omitted the linear part and the constant C which are not critical in our problem. The output 
state will be £{p) oc [G"^^ p] , where the trace is taking on the second part of G^^ and the input state p. The 
completely positive map on the input state will be p^'^ = (£ (8) I)( (V'l )■ The correlation matrix (CM) of the 
Schmidt purification |V') is i;4, 



where 7 is the CM of input state p, (3 



J, 



I are purely off-diagonal, where 



Jn — J 



k=l 



-1 
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Every operators A e B{H) is completely determined by its characteristic function xa{x) := Tr[ylM^(x)] 
[7]. It follows that A may be written in terms of xa as|H] A = n^"^ J^„^ dxxA{x)W{—x). Thus|V') (?A| — 
7r~^" /r<1" d,xXii>{x)W{—x), with Xipi^) — exp[— |a;"^7^x], we assume the first moments of the input state p be 
zero. Hence p'R<3'=7r-4" /jj,4„ dxidx2Xi,{xi, X2)W{-X2)Tr2[G^^Wi-xi)]The CM of p^'?' will be 



ri-ri2(r2 + 7)"'rf2 ri2(r2 + 7)"V 



/3^(r2 + 7)-^r' 



iT-T 
12 



7-/?^(r2 + 7)-V 



(6) 



Where we have denoted 



F = 



Pi 

J- 12 



ri2 
r2 
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and r = (I®A)r(I® A), with A = diag{l, —1, 1, — 1, • • • ,1, —1) is a diagonal matrix which represents the transpo- 
sition in phase space {Xj 



-Pj).The CM of the out output state is 7' = Ti -ri2(r2 +7)"ir^2- 



Now we turn to the trace-preserving Gaussian CP mapslH] which describe aU actions that can be performed on 
p by first adding ancillary systems in Gaussian states, then performing unitary Gaussian transformations on 
the whole system, and finally discarding the ancillas. On the level of CMs these operations were shown to be 
described by 7 i— > 7' = M'^jM + N. The Gaussian operator that corresponds to this operation has the CM|S] 



lim 

r — >C30 



M^ArM + N M^Cr 
CrM Ar 



where Ar — coshrl and Cr = sinhrA. By taking the limitation oi r ^ 00 one has the CM of p^-^ to be 



r' = 



r' r' 

J- 1 J- 12 

^12 ^2 



M^7M + N M'^fi 
P^M 7 



(7) 



This is the final result for a trace-preserving Gaussian map on an input Gaussian state. The reason that we 
restrict ourself to the trace-preserving Gaussian CP maps is from the physical consideration. The result state 
of trace-preserving Gaussian CP map is a state with its CM 7 remains intact in the reference system R (see 
Eq.Q), hence we can compare the output of the Q system with the input state which is keep intact in R 
system. While a general Gaussian CP map will not only change the Q system but also the reference system R 
(see Eq.®). The distortion is some kind of difference between the output and the input. If the input state can 
not keep, the definition of the distortion will lost its basis. Hence we can only define distortion on the basis of 
trace-preserving Gaussian CP maps with a clearly physical meaning. 

One mode Gaussian state: The positivity of p^'^ can be written as the uncertainty relation F' — iJi2,n > 0. 
Due to our selection of the purely off-diagonal /3, we have Ji2,n = Jn® {—Jn)- For a one mode Gaussian state 
input, r' is a 4 X 4 matrix. The uncertainty relation requires det(r' — i J12) > which can be expressed as|10| 

det r; det Pa - det P'l - det T'^ + (1 -f det T[^f - Tr{JT'^JT'^^JT'^JV'^^) > 0. 



The sign before det T'12 now is positive due to J12 = J ® (- 
inequality will reduced to det N' — (1 — > 0, that is 



J). Denote K = detM, and N = M'^N'M, the 



detA^- (1 - Kf > 0, 



(8) 



where we have used det{A + B) = dot A + det B - Tr{J AJ B), M^JM = K and det 7^, = 1. One the other hand, 
the positivity of the output state reads T'^-iJ > 0, thus M^-fM+N-iJ = {^-iJ)M+N-iJ{l-K) > 0. For 
any input state we have j — iJ > 0, the equality can be achieved by pure state, hence we have N — iJ{l — K) > 
which will also lead to Ineq.®. 

The physical meaning of Gaussian trace-preserving map indicated by M and A^ is that N is the additive 
noise and M is a symplectic transformation (rotation and squeezing) and a successively amplitude damping 
or amplification. Let us consider the amplitude damping (as well as amplification) of the channel, which 
is described by det (A/) = K. The amplitude of the signal is damped by a factor of fc = what will 

we do to retrieve the input, clearly we will amplify it back. Or we will reduce the input state with the 
same factor to compare with the output. In these two cases, the distortion operators will be modified to 
d{A,B) = ^[{XA/k - Xb)^ + {PA/k + kPe)^] and d{A,B) = ^[{Xa - kXs)^ + {Pa + kPs)^] respectively. In 
both these situations, if we take all the steps as a whole channel, then we have K = 1. Hence in the following 
we just need to consider the channel of symplectic transformation and additive noise. 

Let us consider the coherent information, which is determined by the symplectic eigenvalues of P' and T[ . 
Now detM = 1, hence (M ® I) are symplectic transformations, which preserve the symplectic eigenvalues so 
that the coherent information. P' can be written as P' = (M^ I)P"(M I), with P'/ = 7 + (Af^)"^ 7VAf^\ 
= P, P2 = 7 correspondingly. 

The problem is to search a M such that the average distortion d is minimized. We have 



F" 
'- 12 



d = ^[TrN + TrM'^jM + Trj - 2TrMf3A]. 



(9) 
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The minimization of d will involve algebra equation of power 4. Let us firstly consider the input of the thermal 
state whose CM is 7sl, with 7s — 2Ns + 1 is the symlectic eigenvalue of the CM 7 and Ng is the average photon 
number of the state. Denote jTrN — Nn, we have 



M = 



\/l + s2 - 52 s _|_ ^ 



VI + S2 - J2 



with s = ^72 — 1/ (27s) = sinliTs, and the distortion operator 

7(— +37s) + iV„. 
4 7s 

The minimal A'^„ is 0, hence the minimal dmin = ji^" 37s), thus we define a canonical distortion dc instead 
of d, dc = d- rfmin, 

dc = (10) 

The coherent information now is determined by the symplectic eigenvalues of T" and F" . The symplectic 
eigenvalues are functions of the trace and determinant of the noise term {M'^) NM~^. Denote S = jdetN, 
T = iTr[(M^)"^ NM-^], the coherent information of the state with CM T" wiU be HHgl 

Ic = gido 9{d, - 1) - 5(d2 - \). (11) 

Where g{x) ~ [x + 1) \og{x + 1) — xlog(x) is the bosonic entropy function, and d^ = x + (Ng + ^)^, di 2 = 
^[x + i ± x^ — 4:Ns{Ns + 1)S], with x = 5 + {Ng + ^)t. The entanglement information rate distortion now 
is the minimization of Ic over all possible noise matrix N with given trace TrN = 4A^„. After the determinant 
and the trace of the noise matrix N are given, we still have the freedom in choosing the off-diagonal elements or 
the difference of the diagonal elements. This freedom and the undetermined parameter k in the matrix M are 
combined into a parameter t (—1 < t < 1) and we can express r as r = 2[Nn cosh(2rs) + t^jN^^^^ sinh(2rs)]. 
Thus 

= m\nlc{x{5,t),5). 
s.t 

The minimization will be achieved when S — iV,^, we prove this by firstly preserving x while increases S. When 
sinh(2rs) > Nn/2, this is always possible by varying t properly to compensate the change of x caused by the 
increase of 5. We have ^^"q^^^'^ = co[/(di — ^) — /(d2 — \)\. where /(a) = 2aTT is a monotonically 

decreasing function and cn = . ^ > 0. Thus ^MT'*^-* < 0, Idx.S) monotonically decreases with S 

increases while preserving x. The minimum is achieved at S = N^, that is, the noise matrix N is proportional 
to the unity matrix. The condition sinh(2rs) > Nn/2 may contain most of the situations. For most of the 
input states (i.e. Ng > 0.012) when Nn = 2sinh(2rs), the values of the coherent information will be at 
S = A^^. We need further to consider the situation of weak signal input states (i.e. A^s < 0.01). There is the 
case that x{0,t) is too small compared with x{Nn,t) for all t. So we need firstly increase x from x(0, ig) to 
some intermediate state with x{Si,ti) = x{Nn,t) while in the x increasing process the coherent information is 
decreased. Let xi,2 = ^[2;^ ± 4:Ng{Ng + 1)6], we increases x and d while keeping X2 invariant, then — " 1' ^ ^ — 
^[/(do — ^) — ^(/('^2 ~ ^) + /(^2 — ^))]- The function / is not only a monotonically decreasing but also a 
downward convex function. Hence in order to prove ^^"''g^^^^'' < 0, we only need to prove do > ^(di +d2) which 
is confirmed by the fact that d^ - ^{dj + dj) = {Ng + + \x ~ \ ^ Ng{Ng + 1) + ^a; > 0. This completes 
our proof. The entanglement information rate distortion of thermal state input as a function of the canonical 
distortion Nn will be 

R'{Nn) = max{0,/c((5 = NI,t = 2A^„ cosh(2rs))}. (12) 
where cosh(2rs) = 1 + 2As(Afs + l)/(2As + 1)^. 
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Figure 1: Entanglement information rate distortion, the trace of input CM is Trj = 3,Ns = for pure state, 
the maximal Ng is for thermal state input. 

For the general input 7, we may further transform the CM F" to F'" by symplectic transformation S (B S, 
S diagonalizes 7 by S'^'yS = 7sl, with 7^ is the symplectic eigenvalue. Thus F'" will be a CM with its 
submatrices are F'/' = 7^! + S'^ {M'^)~^ NM-^S, T'{^ = /3,F^" = 7^! correspondingly. M is determined by 
the minimization of d. The above proving that minimization of coherent information is achieved at 6 = 
remains true for general input 7. The only difference is that we now have r = [A^„r2 + t^N^ — d^Q'^ — 4], 
where Q = = S'^^M. We have O = ^TrM'^^M to be the function of trace and determinant 

of the input CM 7. The entanglement information rate distortion will be 



A special case is the pure state input which contains squeezed states (for coherent states and squeezed 

coherent states the distortion operator should be modified). We have = 2 which is the minimal of Vt. Thus 



i?^(A^„) =max{0,5(iV, + iV„)-5(iV,„i)-5(A^,„2)},whereiV,„i,2 = 1{\/{N„ + l)^ + 4iV„iV, - 1±7V„). Because 
A^s = so that R^{Nn) = for pure states. 



Conclusions and Discussions: We proposed the distortion operator which is quadratic of the canonical 
operators. The distortion operator has a good classical correspondence of mean square error. The distortion 
is the trace of the distortion operator on the joint state density operator. It is an extension of the definition 
of classical distortion. For quantum Guassian state source, we proved that the entanglement information rate 
distortion which is a lower bound of the rate distortion is achieved by Gaussian map under the constrain 
of zeroth, first and second moments. In the language of distortion operator, distortion defined with fidelity 
(1 — Fe) corresponds to the distortion operator of I— l^f) (^fj , where j^f) is the purification of the source state. 
The quadratic canonical operator distortion is more convenient than fidelity distortion for Gaussian state. 

By the physical meaning of the distortion, we rule out the non-trace-preserving Gaussian maps and convert 
the amplitude damping or amplification channels to the standard maps which contain a symplectic transforma- 
tion and an additive noise in the language of correlation matrix. For one-mode Gaussian state input, we proved 
that the entanglement information rate distortion is achieved when the additive noise matrix is proportional to 
unity matrix. The canonical distortion is simply the average photon number of the noise. The rate distortion 
for pure state input is zero. 

One of the most important conclusion we can draw is that the rate distortion function is accessible for 
noiseless case. For any one mode Gaussian input states, the entanglement information rate distortion functions 



7V„0)}. 



(13) 
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at the point of zero distortion are {0) = g{Ns), which is the entropy of the source S{p). From Schumacher's 
quantum noiseless coding theorem^Jj we know that i?(0) = S{p). Thus we have the conclusion that i?(0) = 
RHO). 
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